Distinct predictive performance of Rac1 and Cdc42 in cell migration.

We propose a new computation-based approach for elucidating how signaling molecules are decoded in cell migration. In this approach, we performed FRET time-lapse imaging of Rac1 and Cdc42, members of Rho GTPases which are responsible for cell motility, and quantitatively identified the response functions that describe the conversion from the molecular activities to the morphological changes. Based on the identified response functions, we clarified the profiles of how the morphology spatiotemporally changes in response to local and transient activation of Rac1 and Cdc42, and found that Rac1 and Cdc42 activation triggers laterally propagating membrane protrusion. The response functions were also endowed with property of differentiator, which is beneficial for maintaining sensitivity under adaptation to the mean level of input. Using the response function, we could predict the morphological change from molecular activity, and its predictive performance provides a new quantitative measure of how much the Rho GTPases participate in the cell migration. Interestingly, we discovered distinct predictive performance of Rac1 and Cdc42 depending on the migration modes, indicating that Rac1 and Cdc42 contribute to persistent and random migration, respectively. Thus, our proposed predictive approach enabled us to uncover the hidden information processing rules of Rho GTPases in the cell migration

Robustness of MEK-ERK Dynamics and Origins of Cell-to-Cell Variability in MAPK Signaling.

Cellular signaling processes can exhibit pronounced cell-to-cell variability in genetically identical cells. This affects how individual cells respond differentially to the same environmental stimulus. However, the origins of cell-to-cell variability in cellular signaling systems remain poorly understood. Here, we measure the dynamics of phosphorylated MEK and ERK across cell populations and quantify the levels of population heterogeneity over time using high-throughput image cytometry. We use a statistical modeling framework to show that extrinsic noise, particularly that from upstream MEK, is the dominant factor causing cell-to-cell variability in ERK phosphorylation, rather than stochasticity in the phosphorylation/dephosphorylation of ERK. We furthermore show that without extrinsic noise in the core module, variable (including noisy) signals would be faithfully reproduced downstream, but the within-module extrinsic variability distorts these signals and leads to a drastic reduction in the mutual information between incoming signal and ERK activity

Reconstruction of Insulin Signal Flow from Phosphoproteome and Metabolome Data

SummaryCellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data

Reconstruction of Insulin Signal Flow from Phosphoproteome and Metabolome Data

Cellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data.UTokyo Research掲載「細胞内のビッグデータから大規模ネットワークの再構築に成功」URI: http://www.u-tokyo.ac.jp/ja/utokyo-research/research-news/reconstruction-of-molecular-network-from-cellular-big-data/UTokyo Research "Reconstruction of molecular network from cellular big data" URI: http://www.u-tokyo.ac.jp/en/utokyo-research/research-news/reconstruction-of-molecular-network-from-cellular-big-data

画像データの統計信号解析に基づく細胞運動のフィードバック制御機構の解明

京都大学0048新制・課程博士博士(医学)甲第16731号医博第3679号新制||医||992(附属図書館)29406京都大学大学院医学研究科医学専攻(主査)教授 成宮 周, 教授 山田 亮, 教授 楠見 明弘学位規則第4条第1項該当Doctor of Medical ScienceKyoto UniversityDA

A hybrid in silico/in-cell controller for microbial bioprocesses with process-model mismatch

Abstract Bioprocess optimization using mathematical models is prevalent, yet the discrepancy between model predictions and actual processes, known as process-model mismatch (PMM), remains a significant challenge. This study proposes a novel hybrid control system called the hybrid in silico/in-cell controller (HISICC) to address PMM by combining model-based optimization (in silico feedforward controller) with feedback controllers utilizing synthetic genetic circuits integrated into cells (in-cell feedback controller). We demonstrated the efficacy of HISICC using two engineered Escherichia coli strains, TA1415 and TA2445, previously developed for isopropanol (IPA) production. TA1415 contains a metabolic toggle switch (MTS) to manage the competition between cell growth and IPA production for intracellular acetyl-CoA by responding to external input of isopropyl β-d-1-thiogalactopyranoside (IPTG). TA2445, in addition to the MTS, has a genetic circuit that detects cell density to autonomously activate MTS. The combination of TA2445 with an in silico controller exemplifies HISICC implementation. We constructed mathematical models to optimize IPTG input values for both strains based on the two-compartment model and validated these models using experimental data of the IPA production process. Using these models, we evaluated the robustness of HISICC against PMM by comparing IPA yields with two strains in simulations assuming various magnitudes of PMM in cell growth rates. The results indicate that the in-cell feedback controller in TA2445 effectively compensates for PMM by modifying MTS activation timing. In conclusion, the HISICC system presents a promising solution to the PMM problem in bioprocess engineering, paving the way for more efficient and reliable optimization of microbial bioprocesses

Automatic Quantitative Segmentation of Myotubes Reveals Single-cell Dynamics of S6 Kinase Activation

Automatic cell segmentation is a powerful method for quantifying signaling dynamics at single-cell resolution in live cell fluorescence imaging. Segmentation methods for mononuclear and round shape cells have been developed extensively. However, a segmentation method for elongated polynuclear cells, such as differentiated C2C12 myotubes, has yet to be developed. In addition, myotubes are surrounded by undifferentiated reserve cells, making it difficult to identify background regions and subsequent quantification. Here we developed an automatic quantitative segmentation method for myotubes using watershed segmentation of summed binary images and a two-component Gaussian mixture model. We used time-lapse fluorescence images of differentiated C2C12 cells stably expressing Eevee-S6K, a fluorescence resonance energy transfer (FRET) biosensor of S6 kinase (S6K). Summation of binary images enhanced the contrast between myotubes and reserve cells, permitting detection of a myotube and a myotube center. Using a myotube center instead of a nucleus, individual myotubes could be detected automatically by watershed segmentation. In addition, a background correction using the two-component Gaussian mixture model permitted automatic signal intensity quantification in individual myotubes. Thus, we provide an automatic quantitative segmentation method by combining automatic myotube detection and background correction. Furthermore, this method allowed us to quantify S6K activity in individual myotubes, demonstrating that some of the temporal properties of S6K activity such as peak time and half-life of adaptation show different dose-dependent changes of insulin between cell population and individuals

System identification of signaling dependent gene expression with different time-scale data.

Cells decode information of signaling activation at a scale of tens of minutes by downstream gene expression with a scale of hours to days, leading to cell fate decisions such as cell differentiation. However, no system identification method with such different time scales exists. Here we used compressed sensing technology and developed a system identification method using data of different time scales by recovering signals of missing time points. We measured phosphorylation of ERK and CREB, immediate early gene expression products, and mRNAs of decoder genes for neurite elongation in PC12 cell differentiation and performed system identification, revealing the input-output relationships between signaling and gene expression with sensitivity such as graded or switch-like response and with time delay and gain, representing signal transfer efficiency. We predicted and validated the identified system using pharmacological perturbation. Thus, we provide a versatile method for system identification using data with different time scales